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NATIONAL ADVISORY COBDtETTJSE FOR AERONAUTICS 

-■ ^ 

TECHNICAL MBMORANDUH NO* 1165 

DRAG CORRECTIONS IN HIGH-SPEED . WIND TUNNELS^ 

By H. Ludwleg 

SUMMA^ 


In the vicinity of e .body to a wind tunnel the 
displacement effect of the wake, due to the finite 
dimensions of the stream, produces a pi^resaure gradient 
which evokes a change of drag. In Incdmpresclble flow 
this cliange of drag Is so small, to general, that ono 
does not have to take It Into accoiint to wind-tunnel 
measurements; however. In oomprosslble flow It becomes 
considerably larger, so that a correction factor la 
necessal'y for measurod values • Correction factors for 
a closed tunnel and an open Jet with circular cross 
sections are calculated and compared with the drag 
corrections already known for. high-speed tunnels. 


I. INTI^0DUCTI0N 


The present report deals with the effect of the 
finite dimensions of the stream on the drag of models to 
compressible flow. It will be assumed, that the model 
Is In either an Infinitely long open Jet or an Infinitely 
long closod tunnel in which thore Is a perfectly constant 
velocity without the model installation. Therefore, the 
closed wind tunnel must be slightly tanerod to compensate 
for the effect of the InoreasiJig boundary layer. The 
following tunnel corrections will bo applied for those 
conditions : ‘ • 


(1) When the model' is In- an Infinitely long air 
stream the streamlines near the model are diverted and 


-M- "Wider standskorrektur to HochgeschwtodlgkeltCkon^en, 
Zentrals fur wlssenschaftllches Berlchtswesen der Luft- 
fahrtforsohung das Generali uftzeugmels tors (ZWB) Berlln- 
Adlershof - PB 1955, Gdttlngen, April li|., ISm* 
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the velocities Increased; however. If the' model Is In a 
tunnel of finite dimensions,' thd streamlines can not bend 
outward at the walls of the closed tunnel nor the velocity 
Increase at the boundary of an open Jet. Instead, an 
additional flow appears which produces an Increase In 
the velocity of the approaching flo’./ at the model In the 
case of a closed tunnel and a decrease In the case of an 
open Jet, This correction was extended by B. Gothort 
(1) and A. v. Baranoff (2) to compressible flow. It Is 
designated below briefly, as the displacement correction 
factor. 

(2) On account of the drag, a region^df dead' air 
arises beMnd the body by which the flow 1s deflected 
laterally. Because of the boundary conditions, an 
additional flow again arises at the boundary of -the 
stream which, while It does not produce a ■velocity com- 
ponent at the model In the closed tunnel on the other 
hand evokes an additional velocity at an Infinite distance 
upstream from the model. Since the velocity at an . 
Infinite distance upstream from the model serves as a 
reference quantity, a correction factor Is necessary. 

B. Go’thert (1) calculated It for compressible flows and, 
for short. It Is hereafter called the dead- air displacement 
correction factor. This correction vanishes for free flow. 

Aside from these two correction factors, a third one 
la necessary for the follov/lng -reason. -The additional 
velocity evoked by the boundaries of the stream on account 
of the displacement effect of the dead air produces a 
■velocity gradient and, hence, a pressure gradient. Since 
the model Is then In these pressure gradients. It experi- 
ences an additional drag or forward-acting force . The 
correction necessary to compensate for this change of 
drag will be handled In the present report. It Is called 
the dead- air pressure gradient correction factor. 


II. CALCULATION OP THE COHHECTION FACTOR FOR THE DEAD- AIR 
PRESSURE GRADIENT WITH INCOMPRESSIBLE PLOW 


To begin with, the correction factor Is computed for 
' ,the pressure gradient of the dead air for a closed circular 
tunnel, and a circular open Jet -with Incompressible flow. 
Consider a body of revolution located along the axis of 
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a olosad tuzmel. As a result of it a drag the body has 
a wake, whose effect on the, surrounding air is that of a 
bar with' A drag' area' % = i located behind the 

body, or a source at -the position. of the model with the 
yield (1) 

ft ” 2 ^ ^o 

where c,^ is the drag coefficient, Pjj the reference 
surface, and the flow velocity of the model. The 

quantity Pjj, is called the effective wake area, for 
short . 


To define the velocity grad.ient at the position of 
a source which is located on the a::la of the tunnel and 
substituted for the dead air, the potential of the source 
la taken as (x, r, ’9) whex*e x, r, cp are oyli-ndrlcal 

coordinates and the line r = 0 conlncidoa with the axis 
of the tunnel. The additional flow, with the potential 
along the tunnel v/all, la defined so that 



= 0 


wall 


for a flow jZf = + jZfg 




the quantity { - — is then the desired velocity 

^model 

. gradient at the modal . 


The calculation of the velocity gradient is reduced 
to conventional calculations* Multiply the equation 

+ ^2 toy the tunnel diameter D and differentiate 
with respect to x. 

The expressions appearing in this equation are 
regarded as the potentials of a stream. Inasmuch as 
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represents the potential of a source -with a yield Q, 

D la the potential of a' dipole flow of dipole moment, 

Mp = Q D which Is easily conflriued by substitution and 


dlfferen^atlon of the s our oe. potential, Piirther, 
~ r 6 m ~ ^ along the wall of the tunnel, which 


follows Immediately from the vanishing of 


along the 


tTinnel wall. Consequently, D 




Is the potential of 


an additional flow which cancels the component of the 
dipole flow normal to the wall at the tunnel wall. The 

velocity of this additional flow D ■ - la equal to 


the desired velocity gradient 


T-a except for the 

6x2 


factor D, This additional flow was used. In connection 
with the displacement correction factor already mentioned, 
exactly like a dipole, therefore the calculations that 
have been carried out for this can be applied here . 


Let Pq, Pq, Vq, Mq, and qQ r^ yresent the pressure, 

density, velocity, Hach number, and dynamic pressure In 
the undisturbed flow far upstream from the model, and 
p, p, V, M, and q be the corresponding* values for the 
principal flow plus the additional flow. In addition, 

P - Po = P - Po = ^P» 

At the dipole, or model, of dinole moment Kjj = D Q 
the additional velocity is: 

^^^\nodel ~ V ^ 


according to Gc^thert (1), with a factor depending on 

the tunnel and the method of suspension In the tunnel. 

For the example used - a body of revolution in a closed 
clrouler tunnel - = 1.02. According to the foregoing 

(Av) Is equal to D times the velocity gradient 

BLOdel 
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of a sotiroe Q at the aoiiroe. Therefore 




'Hi = T. 


he cause v = Vn + Av 


model 


Substituting QD for Hp and the expression from equation (1) 
for Q gives 

M = Vq 

^djymodel ” 2 d5 


By Bernoulli's theorem 


dv _ "1 dp 


hy substitution 


2 /dv\ 


2 dv _ -1 dp 
v^ dx “ dx 


model 


'model ” 


■^v ow,.i Py 


This Is the pressure gradient that the wake behind a 
body of revolution In a closed tunnel produces at the 
position of the body. 

If a body Is In a pressure gradient It experiences 

dx 

the following additional drag (5): 

AW = -a ^ 
dx 


where Vjj Is the volume of Idae body and a = 1 + — , 
Vg being ths volume of the apparent mass of the. body. 
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Substitution of the value from equation ( 2 ) gives 

“ ’’t % 9 


& 


d 3 


or 




■°w- 




( 3 ) 


If the model has some type of supports (for example, 
struts or tension wires) with reference surface P. and 

A 

drag coefficient , then the appropriate expression 

XX 

for the pressure gradient is; 

Tv CWji F 


LI +‘*’vA ^A 


model 


Por a strut or tension wire running from the axis of 
the tunnel to the wall, then t Is approximately I.13, 

which can be obtained by extrapolation from Gothert'a 
value for ( 1 ) . 


Por the c^- correction factor it follows that; 


ACr, 




a V, 


. (, , "' A % 

K V ^'h J 


(U) 


.0 


Prom equations ( 5 ) and (I4;) it is. evident .that, the cor-, 
rectlon factor for the dead-air pressure gradient Is made 
considerably larger by high drag at the sup.nort. 

The derivation for a free .stream can be carried 
out In exactly the same mazmer. In this .caa’a it is only 
necessary to require that no change of velocity occur 
at the boundary of the stream. Then exactly the same 
equation as (Ij.) Is obtained, the only difference being 
In the value for which Is =- 0 . 265 , for a body 

of revolution In an open Jet. 
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Thorefore, the dead- air preaaure gradient In a 
cloeed tunnel gives a positive value for tjaat Is, 

an increase In 1 irag. In' a' free stream It dausds a 
decrease In 0^. Consequently, the drag values measured 

In a closed tunnel- must be reduced by the correction 
factor for the’ dead- air pressiu^e gradient to obtain the 
values corresponding to the free medium. Conversely, the 
.values measured In the free flow must be Increased. 

III. APPLICATION OP THE RESULTS TO COMPRESSIBLE PLOW 
(a) Definition of the Effective Wake Area 


The effective dead-air surface for Incompressible 
flow was ; 


= I 


The corresponding exuresslon for compressible flow 
will now be defined. According to the momentum theorem, 
the drag of a body Is given by: 

W = pv(vo - v) df 

when the Integration over the wake is carried out In a 
section behind the body where the static prossxire has 
again reached the value p^. Further, it Is assximed 

that this section la located so far behind the body that 
Av = Vq - V « Vq and’Ap = Pq " P « Pq 

Approximately, then 

Po 


or 

2 



( 5 ) 
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For the effective dead-al-r aurfaoe- Pu, and^ therefore, 

for. the siirfaoe over which thd flow is deflected by the 
wake: ■ ^ ' 


= Po ’'0 




Vq - pv) df 


or making the sazne approximations as before 




A3sumJ.ng no heat transfer at the model, it follows 
from the energy theorem 


■ where i represents the heat content. Since pv df = dm 


is always positive, the expression 


i^o* ~S~~ 


must vanish across the wake, on- the average. As an 
approximation, it can be sot equal to zero for each 
streamline, which then means that there is a constant 
stagnation point temperature in the v/ake . This assumption 
is also used in the evaluation of loss of momentum 
measurements (q.). Therefore, 


lo + 


, v2 
= 1 + -^ 


Proceeding from the as simp t ion that p = p^, and 
introducing the quantities Av and Ap while bearing 

in mind that 1 = ~ I ' T « 


^ .(k - 1) 


By substitution in equation (6) 


= 1 


1 + (ft-1) A* 


■1/s 
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and referring to equation (5) 

% = 5 E * t" - “oH 


- (7) 


Therefore, the effective wake area Is greater In 

ooinipresalhle flow for the same . than the- value In 

Inoompreaalble flow by the factor |~1 + (k - L) . 

This la due to the fact that the wake In compressible ' 
flow has a lower density, as well as a lower velocity. 

(b) Application of the Pomrula for the Correction Factor 
for the Dead-Air Pressure Gradient with the 
Aid of the Prandtl Theory 

To apply equations (3) and (4) to oomproaslble flow 
the Prandtl theory will be used In the fona of the 
streamline analogy (5). For the present problem It 
reads as follows: 

Assume a model In flow at a Mach number Mq In a 
tunnel of diameter D; In a comparison tunnel with 

diameter D ^ with a model in incompressible 

flow whose dimensions Including the wake, have been 

reduced by 0. - at right Mngles to the direction of 

flow and remain unaltered in tVio direction of flow. The 
additional longitudinal voloclties appearing In the 

‘ 1 

compressible flow are then x times greater than 

1 - M 

o 

for the corresponding points In the Inoompreaalble 
comparison flow. 

However, at corresponding points in the comparison 
tunnel the additional velooltlee are '■■■■■■ tinea 

Vi - V 

larger on account of the displacement of the dead air 
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than In the original tunnel with Incompressible flow. 


since the model Is . 


- “c 


longer in the Incompressible 


comparison tunnel relative to the diameter, while the 
cross-sectional dimensions of the model and dead air are 
equally large relative to the tunnel diameter. Therefore, 
for the same dead-air displacement in compressible and 
Incompressible flow In the same tunnel, the velocity 


gradient increases by the factor 


a - 


Besides 


since the dead-air displacement for the same values 
of c^ Increases by t + (K - 1) the velocity 


gradient Increases by 


the same value. 


1 + (k - 1)Mq^ 

(i - 1.35/2 


If the 0 retains 
w 


Consequently, equation (2) talces the following fomr* 
for compressible flow . 


2 

V 


model 


_ \ 1 + (K -'l)*^o^. 

- K T' • ■ * 


r.5 (1 - 


( 8 ) 


Beraoulli’s equation gives 


dp - o 
. q “ ” ^ V 

for compressible flows'. It follov/s fi’om this, that. 


1 

<1 


_ V 


‘model 
and when we again set 


d 3 (1 - K35/2 


dp 
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there follows 


or 


AO. 




*M 


_ g V^M * (« - DHq^ 
d3 . H^2J3/2 

(1 * Cst'^u^^ * 

d3 (1 - m„2)5/^ 


( 9 ) 


( 10 ) 


Equation { 9 ) holds for a model if the drag effect of the 
supports la Ignored, while this has been taken Into 
account In aquation (10). The values to be assigned to 
the factor a for compressible flow are discussed In 
the ne^t section. 


(o) Bodies In Compressible Plow with a Pressure Gradient 
(Definition of the Factor a) 

It Is a well known fact that a body of volume 

In an Incompressible potential flow with a pressure - 

gradient ^ experiences a drag 
dx 

W = -a ^ 
dx 

V 3 

According to 0. I, Taylor (3), a = 1 + ~ where Vg 

is the volume of the apparent mass of the body.' For 
slender, streamlined bodies Vg'« Vj^, .consequently, 

a Is approximately equal to 1. The values of the factor 
a for compressible pubsonlc flow with'pressure gradients 
are obtained by application of the energy theorem. 

Consider a, ccanprosslble stream of Infinite extent 
with a pressure gradient In- the direction of flow. The 
pressure.,, density,^ velocity, heat content* and Internal . 
energy of the ^disturbed stream are denoted by ff, p, 

V, T, and u, respectively. Vi/hen a body la In the stream, 
let p, p, y, 1, and u represent the disturbed values in 
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th6 vicinity of the body. Purther, aaatimo that there la 
perfect flow around the .body. Then, at greater dlatancea 
from the body p = p; p = ^, etc. To continue, the 
Bernoulli equation 


1 



C-onatant 


holda for all points In the flov; in a reference sj^teri 
fixed In the body. How Imagine a control surface around 
the body with a large distance uet^Tben the two, displace 
the body by an amount “(Ut ' ^wi'Sh the control surface 
rigidly fixed and apply the -energy the of 5k to the -flow. 


The energy removed from the flow by the diaplacement 
la then equal to «7 dx, where V.' is the. drag on the body. 


P ^ ev 1 


Is the energy of the gas 


enclosed in the control surface. If Is tho volume 

enclosed by the control surface and if p = 0 at the'' 
positions of the model, in evaluatlJig tho integral, 

Tho change in the enclosed energy produced by the 
displacement of the body la then; 



The increase in energy with^ the control space dus 
to flow thi'ough tho control surface curing t)i^ ^IsnJ.acement 
is 



when ti is an Instant prior to the st^rt- of thedls- 
place’ment and .tg^, an instant after the displacement 

at which the stationary state .is ag.",ln reached , P., is 

. • , ■ ■ . 1 ■ 1 1 ■ 

the "control surface- and Vjj 3.3 th^ component of the 
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velocity normal to the control surface, Dvu?lng the . . 

.dlsplaoemejQt .of '.the. mods! the JBemoulll equation does’ 

not anply In the form: ‘1 + -rr = T +• -r = Constant since 

2 -2 , . 

thib flow In the present reference system (at rest> referred 
to the ■original position of the hody) Is no longer station- 
ary, because of the movement of the model. -A closer- 
examination shows, however, that these devlatliOns- fall 
off so rapidly with Increasing distance- from the modal 
that for very large dlsttoces -of^the control si:^face froto 
the model, nevertheless. It Is possible to wrlt'e 

l-i-~=l+~ = Constant Ahd then .bring It out In 
front of the Integral. Therefore 

V -2V fff 

Jt^ pv^ Jti yj p^n dp- dt . 


The integral on the right side simply represents 
the excess of Incoming over outgoing mass between the 
times t^ and TherefdT*e It foll^-ws from the equation 



According to the energy theorem the energy of the 
stream- taken up durlhg the displacement must equal the 
loss of energy of the gas enclosed In the control space 
plus the excess of energy entering the control surfaces 
over that leaving. Therefore 1 



Since. the undisturbed quantities p, p, v, T, and u are not 
liffected by the displacement 
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(= . 


1^=0 


‘^LTk ^ ‘IVJ = 0 and p =0 

whore p, u, and v at the poaltiona of the model are 
ae signed the. values that prevail '.when no model la present. 
Prom this If follows- _ 


W dx = •» d 


TJ _p (" * f). -p (" * #)] 
* ^ ^ [if/ (p - p) '*’'1 


' J 

The volume Integrals which covor the entJ.re volu’vo "Vk 

enclosed by the control surface are split into two 
Integrals^ pne of which .extends ovor Vjj the space 

occupied by the model, while the other extends over the 
apace Vg; - Vjj occupied by the fluid while observing 
Immediately that by assumption p = C in the space Vjj. 
Therefore _ ^ 


W dx = d 


//p (-f) 

TW ' ' 


L vm - 


fii 

Vk - Vji 


(^*t) *(|f/p' 
- p("^t)] 


O u ■ 
Vk • Ym 


(p - p) 


(11) 


Essentially, the only contributions to the value of both 
integrals ever dome from the immediate vicinity 

of the model since, p-^p, u— ^u, and v — >v at points 
more distant from the model. Therefore, for sufficiently 


large values of 


the Integrals become Independent of V_, 


/ - 
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Assmlng only a small presaure gradient In the 
xmdisttirhed flow. It oan then further he assumed as a 
good approximation that the discrepanoles between undis- 
turbed and disturbed flow in' the case of flow with a 
pressure gradient are equal to the dlaorepanoies betvfoen 
the -same quanti^ties in a parallel flow with the same Mach 
'number. Now the last two integrals of equation .(11) can 
be evaluated for corresponding parallel flow. 


• V/lth Po, Po» Vq, lo, and Uq denoting tho undisturbed 
magnitudes of a parallel flow 



(P 



- p) dV 




( 12 ) 


in which the left side is to be evaluated for flow with 
pressure gradient and the right side for a. parallel flow, 
and Po, Uq, Vq must be taken equal to the values of p, u, • 
at the location of the model. 

To. define the two integrals on the right side one 
proceeds as follows • On accelerating a model in incom- 
pressible. flow there is a! drag W where 


W = p 


dv 

dt 


* __ u 

in which ~ is the acceleration with respect, to the 
dt •• 

fluid Vg the volume of the so-called apparent mass. 
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Correspondingly, for compressible flows 

■ *» = Po Vg f (Ko) g (15) 

where Is the undisturbed density at a distance from 

the model, V the value of the apparent mass for Incom- 

Q 

presslble flows, — the acce 7 ueratloii relative to the 

dt 

fluid and f (M_) a function of the Mach nvunber. The 
magnitude of the function f .will be given closer 

attention later on. 

Consider, now that the model Is In a parallel flow 
and moves along with the .stream. Then, p = p^, p = p^, etc, 

.In the entire apace occupied by the stream. The flow Is 
viewed from a reference system In which the undisturbed 
velocity la v^. At a great distance from the model there 

is a perpendicular control surface. Now^ let the model 
be slowed down to rest, gradually €ind, again, apply the 
energy theorem to the fluid In the control space . 

During the process of slowing down the following 
energy la taken from the flow 



where 


and t. represent the place and time 


1 * 2 ' 1 ' “2 
of the beginning and end of the displacement. 


Is 


the relative velocity of the model with respect to the 
undisturbed flow. Prom this one obtains, thon, by 

transf brmaljlon and tekUng Into account that 7— = M 

Vo o 

and — = 
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= Vs Po 


f(M) (Mq - M) dM (llj.) 


which by setting 


f. (M) (1^ - M) dM = g(M ) 


"O do 


becomes 


( 15 ) 


W dx = Vg Po -f- g(Mo) 

2 

The change In the energy enclosed In the control apace Is : 

/II [" ^ 

R. I'i t 

where the integral is taken ov;':"' the vicinity of the 
model and Vj- is chosen large enough so that further 

increases do not affect the integral. 

For the excess of energy added over that taken away 
through .the control surface during the braking, one has: 




dt dF 




(p - p„) av 


% - Yh 
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v2 

and here, too, i + — approaches 1^ + so rapidly 

with Increasing separation frQm the model that one can 
write this expression In front of the Integral, for a 
stiff Iclently large control surface. 

t 

The energy theorem then reads: 

V 2 P 

'^S Po -f- 8 («o> = - 



^ + f) - Po + -f-) 


rT/ 


Vk - \ 


V (p - Po) dV 


However, the right side Is exactly equal to the right side 
of equation (12). Prom this throuch auhatltutlon in 
equation (11) the following Is obtained: 


W dx = d 




'M 



+ d Vg p Y s 


(16) 


d£ 


Assume now, that the pressi\,.'o gradient t*- la small 

“X do 

and constant along the chord of the model, then 
dv 

—etc. must also be conaioerea constant along the model 
dx 
and 


W dx 


= h P (" •" ?) ^ to 

+ ^ [p T s . 


dx 
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OP taking Into account 

. n 

‘ ■ Y-.- ^-'4; & J. + *— ■=■ Gonstafit- -- 

P = • ■ 2 . ■ ■• ■ 


' W 


= - to % + S P T8.<®>_ 


( 17 ) 


For adiabatic flow one has 

d 


(f- t) 


dx 


( - f)i 


^ - M ^ 

dx " “ 2 V 


K - 1 ~ 


ijja) -^35 

/ _ ^ dx 

^ n *T 


P 2 


w - - ^ V - 

” dx - 


^ Vo 
dx ''3 


By substituting In (I 7 ) 

g (H) + I (1 + S' (H) 

and, to get a :nore' comnaot expression. Introducing 

h m = (1 - t)s I s* 


then 


which gives- • 


” dx 


(vm 


a = 1 + 


for the factor a In equations (9) and (10'). 


(E) Vgj 

( 19 ) 

1 

(H) 

(20) 

} and (10'). 


definition 

of the drag 
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and a pressure gradient If the acceleration .drag 

of the same model fop the, range of Mach numbers between 
-9 and S. and, therefore' the • function f(^) In 
equation ( 13 ) la known. 

For Incompreaslblo flow (M = 0), f(0) In equation (15) 

become a 1, From that one likewise obtains 1 for the 
value of h(0), and e.quatlon (I 9 ) goes over to the 
familiar formula developed by Taylor 

I 

The additional drag experienced by a wing which 
satisfies the assumptions of Frandtl's theory when It Is 
In a flow where a pressure gradient exists will be 
Investigated. The first thing to determine la the value 
of acceleration drag for such a wing that will give the 

function f(H) from which h(M) can than be obtained, 

• * * * 

To’ begin with, let the wing liave a stationary motion 
with a velocity v^ In a medlian that has the undisturbed 

values' pQ, pjj, 1^ and Uq at a distance from the model. 

In a reference system with axis fixed In the wing the 
potential of the stationary motion then bocoinos: 

“ ^o * 7 , z, M) 

where qjg approaches sero far from the model. Now If 

the wing moves with the same velocity and a small acceler- 
ation b, then the velocity field can be considered quasl- 
statlonary. The assumption that the acceleration b be 
small Is necessary here in compressible flow In constrast 
to the incompressible flov/, since all disturbances must 
spread out with a finite velocity and since at larger 
accelerations they produce the larger pressure differences 
appearing between accelerated and unaccolerated motion, 
as well as density differences and the corresponding 
velocity differences. Then the appropriate potential ^ 
la: 
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where Vq l3 the velocity attained exactly at the tine 

t = 0. How to find out how thd 'difference of pressure 
appearing on the wing surface, at the time t = 0, 
between the stationary and accelerated motion depends 
on the Mach number; for this pressure difference certainly 
yields the acceleration drag. 

For the stationary motion with the assumption of 
adiabatic flow, Bernoulli's equation reads: 


1 


S 





( 21 ) 


Applying the more general form of Bernoulli's 
equation to the motion that Is not stationary: 




v2 


^ + P- TT = Constant 

t 2 


The heat content 1 can be substituted for the 
pressure function P again. TT Is the potential of a 
force field that might bs present. Since the reference 
system fixed In the body is an accelerated reference 
system, it is necessary to set U = bx. For the Instant 
t = 0 the general Bernoulli equation tlien reads 

^ <Pq (x, y, z, M) + 7 ^ ~ + ^+ i-bx = bx + 


du 


+ 1 - bx 

o 


or 


^o ’S 


dep, 

k 

dji a© 2 


r- f a (*. 7 , K) + ^ ^ + 1 = -^ + (2^) 


2 


Since the flow is quaalstatlonary with respect to the 
velocity field 

’ 3 ^ _ t2 
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From this and equation (21) and (22) It follows that: 

la - 1 = A1 = ”S (*. y* ^ ^ 

For adlahatlc flow ^ = dl Is valid. Since a small 

P 

acceleration and slender model have been assumed, as an 
approximation, therefore one can set: 



A1 


Then from this it follows that: 



The quantity <Pg Is the potential of the Increased 

velocity appearing at the profile, Thei’efore, according 
to the Prandtl theory the following holds at the profile 
surface 


‘Ps 




incompressible 



(1 - m2) 



Incompressible 


Consequently, Ap becomes, at the profile surface 


A -• ^ 

^ ^ ^o "^o ^incomoressible 


S o * ■ ■ (p 

Vq inconroressible 


(1 - m2) 
1 


1/2 


m2 


(1 - 


M^) ■ 


(1 - ll2y 


5/5 


The acceleration drag that appears Increases by 
^ • Therefore, according to equation (13) the 
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following holds for TCJ/Iq): f (M) = ^ 

Prom that and equations (15) and (18) the functions 
g(M) and h(M) also follow 


g(M) 

= - r 

L- 

L - (1 - m 2 )^/^^ 

(23) 

h(M) = (1 - 

^ g(M) 

+ M.(l + « - g«(M) 

(2h) 

\ 

2 V 2 / 


The function 

h(M) 

is shown in figure 2 . 

One sees 


that the factor a = 1 + h(M) — — appearing In equation (9) 

and ( 10 ) is Increased further by the effect of compressl- 

Vs 

bllity. For slender bodies for which is small 

compared to 1, one can therefore also set a = 1 for 
compressible flow as an approximation. 

The argviments carried out to define the potential 
drag of a body in compressible flow with a pressure 
gradient assume that the- flow is a perfect potential 
flow. The validity of the Prandtl theory, was only 
assumed for the definition of the acceleration drag to 
learn the function f(M), 

Even after the validity of the Prandtl rule has been 
exceeded, with the appearance of local supersonic fields, 
however, the flow can still be considered as loss-free 
potential flow, if the subsequent passage of the flow 
into the subsonic region takes place adiabatlcally and 
also if, as a coarse approximation, the flow remains 
loss free in the face of weak shocks (naturally, neglecting 
the wake due to friction and separation of flow). The 
increased velocities appearing at the profile circumference 
are underestimated, in this case, when obtained by Prandtl’ 
theory. Therefore, by the same arguments advanced in the 
derivation of the function h(M), in the application of 
equations (25) and ( 2 l|.) at Mach numbers exceeding the 
critical value, the Correction factor is somewhat 
unde re s t imated . 

For slender bodies of rotation in incompressible 
flow Vs is very small compared to Vj/j# so one can 
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set a = 1 as a good appr oximat 1 on . Now, since the 
Induced velocities Increase considerably slower for 
the body of rotation than for the wing, the function 
h(M) deviates from 1 . still less, so tliat one can 
set a = 1 as a good approximation In compressible 
flow for bodies of rotation. 


IV. LIMITATIONS- IN TEE APPLICATION 
OF THE COHRECTIOIT TACTOR 


So far, the dead air has s3.mply been replaced by 
a source, at the center of gravity of the model, df such 
a strength that the displacement it produced downstream 
from the model was equal to the c.ead-air displacement. 

The velocity gradient of the additional flow was then 
defined at this source on the assumption of constancy 
over the entire chord of the model. Actually, the 
substitute source should be split up Into many small 
sotxroes and suitably distributed to fit in along the 
model aiid e-ventually the foirward part of the wake , 

Tlien the additional flow of all these . sources should 
be determined and from this by integration the velocity 
gradient w 2 ilch could vary along tlje , model. 

To check to what extent .the above- slmpllflcatlcns 
are permissible, consider figure 5 . 

Here the variation of the.. additional velocity along 
the tunnel axis is given for a source Q on the tTumel 
axis in' incompressible flow. According to figure 3 the 
curve for the closed tunnel is -defined for a dipole . 
from the additional flow by Integration. The values for 
the additional flow for a dipole are taken from a calcu- 
lation l5y V, Baranoff (2) for- the displacement correction 
factor.' Suoh a single source produces an approximately 
cohstant velocity gradient of the additional flow for 
approximately 0.7 tunnel radius, both upstream and 
downstream. Therefore, if. 'the model .is .short (less than 
0.7 tUniiei radius ) it • Is iinpateriai hpw. the single 
sources- ore distributed iij! 3 .iijj 3 .i,th 6 model- because the 
resulting velocity grAdlent la . alwa;^ saiM as fcrr a- 
coBiblned source of equal total strength. However, since' 
the principal part of the voliane of a longer streamlined 
body Is in the' center of .the model and, in addition, the 
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substitute sources are to be applied to the central 
-pqrtlon of ^the ^body, essentially,, t^ formula derl^d 
for Incompressible flow can bo used! up 'to' mddo'l lengths 
of one tuzmel diameter. 


The variation of the additional: velocity for 
compressible flow for the ^ame breadth of wa^ Is obtained 
by the Prandtl theorem by Increasing the otdlnates by a 


factor of 



and decreasing the abscissae by a 


factor of 



Therefore, the model length for which calculations 


can be made with constant velocity gradients is 




shorter for compressible flow than for Incompressible 
flow. If the model exceeds the permissible length, the 
amount by which the correction factor from the derived 
formula is too large can be estimated from figure 3* 


The variation of the additional flow of a source 
in a free Jet with incompressible flow also appears in 
figure 3 • The curve is taken frcaa the work of 
D. RCfchemann and P. Vandry (6). One obtains the variation 
for compressible flow by the same distortions as for 
the closed tunnel. It is apparent, that the correction 
factor for tho free flow has tho opposite sign from the 
closed tunnel. In addition, the correction factor is 
only about a fourth as large as for the closed tunnel. 

The length over which tho pressure gradient can be 
considered constant is likewise smaller for free flow. 

By means of figure 3 an estimate can be made of the 
amovuit by which tho correction factor formulas indicate 
the correction factor too large for longer models. 

As stated already, the correction formula for 
free flow with the same size models In tiuanels of the 
same size Is approximately one fourth that for the 
closed tunnel. With long models, for which the validity 
of the correction formula has already been exceeded, 
since It can not be evaluated over the- entire length of 
the model with constant pressure gradients, the effective 
correction factor with free Jet Is considerably smaller 
than- a fourth of the correction factor for a closed tunne^ 


z6 
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^or the fomtulas which apply for constant pressure 
gradients lose their validity with a^ free jet even with 
smaller models, and Indicate the corrections too large. 

V. OOrjPARISON OF THE CORRECTION FACTOR DERIVED WITH OTHER 

drag correction factors at high hach numbers 


According to Qothert (1) the displacement correction 
factor furnishes the following values for the dynamic 
pressure, or Oyf, without taking into account the support: 



On comparing these cori’ection factors with the correction 
factor for the dead-air pressure gradient (equation (9))» it Is 
seen that for Incompressible flow (li = 0) the 
correction factor for the dead-air pressure gradient up 
to the factors a or which are nearly 1 for 

slender bodies, is half as large as the displacement 
correction factor. Moreover, for Incomprosslble flow 

Vs 

a = 1 + — and Is a factor with which the substl- 

Vm 

tute dipole strength can be obtained from the volume of 
a body and the velocity of flow near It. As Qlauert (7) 
who defines somev/hat differently has shown, a = 

exactly for Incompressible flow. 

The Increase of the correction factor v/lth llaoh 
number is different, however, for the two correction 
factors. 

In figure 1 the ratio of compressible to Incoinprosslble 
correction factors for conditions alike In other respects 
1s plotted against Mach number. The correction factor 
for the dead-air pressure gradient rises considerably 
faster than the displacement correction factor, so that with 
M S 0.35 tho dead-air pres sure -gradlont correction factor 
la already just as Isrgo as the displacement correction 
factor, T/hile at M = 0 It was only half as large. 
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According to'-Qothert, the dead -air displaoeraent 
correction -factor furnlahea' for a modal without taking 
Into accottait the aupport 

Aq _ ^ 

ir (1 - MqS) 

or taking In-tO account the apreadlng out of the dead 
air thrcugh comprea3lhy.lty end Ill(a) 

^ V - y) [1 - 11 «o^ 

Si n d2 (1 - Mq2) 



The rise of this correction factor .vith Mach number 
la likewise shown in figure 1, It rises even slower with 
Mach number than the displacement correction factor. 

Since the dead-air displacement correction factor for 
incompressible flow and slender bodies with small c^ 

la smaller in general than tlio other two correction 
factors, it lags behind these even mo7?e at higher Mach 
numborc. 


Through the drag of the supports, the do ad -air 
pres sure -gradient correction factor can be increased 
considerably as shown in aquation (10). It can then 
considerably exceed the displacement correction factor 
if the s-upports have very large drag for relatively 
small volume (for example tension v/ires or struts, on 
exceeding the critical hach number). 


As an extreme example, consider a streamlined body 
of rotation with a diameter of 100 millimeters, and a 
volume 2 liters, in a closed tuniiel with a diameter of 
1 meter, Tlie body is suspended oh six strong l-mllllmeter 
gauge sonslon wires which nm perpendicular to the flow 
from the model to the wall. The drag boeffiolent for 
.the model is = 0,05, "the tension wires it is 

Cw = 1«3 and the Mach number is 0.3. Then oqviation (10) 


gives : 


Ac 


WK 


'm 


= 0.13 


2a 
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This means that, in this case, the drag. Is In error hy 
15 i^rcent .because -of the dead -air pressure gradients. 
Displacements correction .factor and. dead**alr displacement 
correction factor only give 1.5 percent and 0,65 percent, 
respectively, according to Gb'thert. 

This example shows tha^t attention must be given to 
this problem In high-speed tunnels In the measurement of 
models with low drag, to keep the drag of the supports 
very low sc that dead-air pressure gradients will not ■ 
give false drag indications. Tension wires. and, struts 
too, for which the critical iriach number la reached 
rather early are, therefore, extremely disadvantageous. 

• VI. SmiARY 


It has been shown that in addition to the usual 
drag and Kach number correct.lon factors for high-speed 
ttmnels which take into account ciiange of velocity at 
the position of the model due to the model and dead- 
air displacement, a further drag correction Is necessary. 
This correction facton designated the dead-air prosstire- 
gradlonb correction factor. Is based on the fact that 
a pressure gradient arises when the dead air departs 
from the position of the model, which produces an 
additional drag or forward-acting force . It he-s been 
shown that this correction factor is of the same order . 
of magnitude as the other two correction factors and 
that It can even be considerably larger in special cases. 
Ii^ addition, it rises the fastest of all correction 
factors with an Increase in Mach number under conditions 
that are otherwise alike. 


The following formula has been derived for tlie dead- 
air pressure gradient correction in the closed circular 
tunnel and the circular open Jet; 




T 

V M 


'M 


, . “''a 

1 + ■' -■ ■■ 
T c_ p 




1 + (K - 1) M 


2 

o_ 




d5 (1 - Mq2) 


5/2 


In this, a is a factor which depends on the sliape of 
the model and 1s somewhat larger than 1 for streamlined 
bodies. For streamlined bodies of rotation a = 1 Is 
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ft' good 'appPoxlmatloxi. For a wing 


a = 1 + IKMq) ^ 

■ V|j 


wbsre- 




is the. yoltuna, of ^h^‘ apparent maa'a oT the 
wing In inoomprehalhle' flow and the volume of’ the 

wing; ■ h(Ho) la- a function of the Mach huniber which la 
ahown In figure 2 . - ' 


and are factors which depend on the 

tunnel shape and the Installation of the model, that 
Is the mounting In the tunnel. For "bodies at. the center 
of the tunnel with small cross extent relative to the 
tunnelr diameter, = 1,02 In a closed circular tunnel 

and T.^ = -0,263 In a circular open Jet under the 

' same conditions . 


According to .GSthert In the closed tunnel 

for -a wing whose wing span 3 Is no lbngar> small in 
comparison with the tunnel diameter D the factor t.^ 

become 3 


B/D 

0 

0.25 

0.5 

0,75 

^v 

1.02 

1,014. 

1.06 

'1.10 


By extrapolation to B/D = 1 , 2 l'.l 6 Is obtained 

for suspension devices (struts aiid tension w.lres) 
reaching from the center of a closed tunnel to the tunnel 
boundary. The quai^tlty Is the volume of the model 


and Owfi» Pm» and F^ represent the drag coefficients 


M 


A 


and, respectively, the re'ferenoe surfaces of the model 
(suDscrlpt M) and the support (subscript A). 


The function 


1 + (K T l)-^o 


(1 - Mq2) 


5 A 


Indicates how many 


times -larger the dead-air presstiro gradient is for a 


Mach number 
fig. 1 ). 


Mq than In Incompressible flow (compeo’e 
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In conclualon. In determining drags througli loss of 
mo&^ntumi measurements, the .dead-air pres sure -gradient 
cdrrebtlon factor disappears since the additional drag 
due. to the pressure gradient of the wake is a p\ire 
potential drag which is not associated with any loss 
of flow ahd, oh that account, is not noticeable in 
the wake. 


Translated by Dave Felngold 
National Advisory Committee 
for Aeronautics 
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Fig. 2 
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Figure 2. h (Mq) for a wing. 
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Fig. 3 



Figure 3. Additional velocity for a source (Q.) along the tunnel axis. 
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